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Abstract 

We employ an effective-field theory with correlations in order to study the phase 
diagram and ground-state magnetizations of a selectively diluted Ising antiferro- 
magnet on triangular and honeycomb lattices. Dilution of different sublattices with 
generally unequal probabilities results in a rather intricate phase diagram in the 
sublattice dilution parameters space. In the case of the frustrated triangular lattice 
antiferromagnet the selective dilution affects the degree of frustration which can 
lead to some peculiar phenomena, such as reentrant behavior of long-range order or 
unsaturated sublattice magnetizations at zero temperature. The selectively diluted 
Ising antiferromagnet on the honeycomb lattice is obtained as a special case when 
one sublattice of the triangular lattice is completely removed by dilution. 
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1 Introduction 



It is known that the pure Ising model with nearest-neighbor antiferromagnetic in- 
teractions on a triangular lattice shows no long-range ordering due to a high degree 
of frustration present [HE]- However, the typical high degeneracy of the ground 
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state can be lifted by a small perturbation, such as an external magnetic field or 
the presence of quenched vacancies. The former case has been shown for a certain 
range of the field values to result in a phase transition between the ferrimagnetic 
phase with two sublattices aligned parallel and one antiparallel to the field (4-Tt) 
at lower temperatures and the paramagnetic phase in which all spins are aligned 
parallel to the field (ttt) & t higher temperatures [3H5]. A characteristic feature of 
the ferrimagnetic phase is a broad m = 1/3 magnetization plateau. On the other 
hand, introduction of quenched magnetic vacancies in zero field locally relieves frus- 
tration and supposedly leads to a spin-glass order [BHH] . A recent Monte Carlo study 
of magnetization processes in the system with both the field and uniform magnetic 
dilution showed that, for example, even a small amount of vacancies can deform 
the broad frustration-induced 1/3 magnetization plateau into a stepwise curve [DJ. 
In our recent study of the critical behavior of such a system we found that the 
interplay between the applied field and the frustration-relieving dilution results in 
peculiar phase diagrams in the temperature-field-dilution parameter space, involving 
multiple reentrants [ID] . Even more interesting is the case when, instead of uniform 
diluting, nonmagnetic vacancies are introduced into the lattice selectively. Namely, 
if the lattice is split into several (in the case of the triangular lattice three) sub- 
lattices which are populated by the vacancies with different probabilities, then the 
geometrical frustration can also be relieved globally, giving rise to long-range mag- 
netic ordering phenomena even in the absence of the field. Kaya and Berker [TTJ have 
shown that if only one out of three sublattices of the triangular lattice is randomly 
diluted then a long-range order can develop in the remaining two sublattices already 
at relatively low concentrations of the vacancies, while the diluted sublattice remains 
disordered at any concentration of the vacancies. Also Monte Carlo simulations [T2] 
have supported this finding. 

The goal of the present study is to extend the study of Kaya and Berker [TTJ in 
order to establish a three-dimensional phase diagram of the system with two sub- 
lattices jointly diluted but with different probabilities. Since the concentrations of 
the magnetic sites of the diluted sublattices represent two independent parameters 
that are able to control the degree of the frustration of the entire system, we can 
expect a nontrivial behavior displaying qualitative differences in different regions of 
the parameters space. Moreover, if one sublattice is completely diluted, i.e., effec- 
tively removed, the lattice reduces to the honeycomb type. Thus, by diluting the 
second sublattice we can obtain the phase diagram of the selectively diluted Ising 
antiferromagnet on the honeycomb lattice as a special case. 



2 Formalism 

Let us consider the Ising model on a triangular lattice, described by the Hamiltonian 

^ = -JE^M> (!) 
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where Si = ±1 are the Ising spin variables, J < is the exchange interaction cou- 
pling, and is the sum extending over all nearest neighbor (NN) pairs. are 
quenched, uncorrected random variables chosen to be equal to 1 with probability p 
when the site % is occupied by a magnetic atom and with probability 1 — p other- 
wise. Then p represents the mean concentration of magnetic atoms. 

A uniformly diluted system, in which the vacancies are evenly distributed over 
the entire lattice, was investigated in the presence of an external field in our recent 
study [10] by the use of an effective field theory (EFT) with correlations (for review 
see, e.g., [IS])- In order to include the geometrical frustration effects within EFT, 
the triangular lattice was decomposed into three interpenetrating sublattices A, B 
and C (see Fig. UJ, such a way that spins on one sublattice can only interact with 
their NNs on the other two sublattices. Thus all the NN interactions are accounted 
for and the frustration results from the effort to simultaneously satisfy all the mu- 
tual antiferromagnetic intersublattice couplings. It should be noted here that the 
EFT with correlations is based on the differential operator technique introduced 
into the exact Ising spin identity and is almost as simply formulated as the standard 
mean-field theory. Another advantage of this method is that it enables systematic 
inclusion of the effects of spin correlations. However, in this work we adopt the sim- 
plest approximation which neglects correlations between different spins but it takes 
the single-site kinematic relations exactly into account through the Var der Waerden 
identity. It is worthwhile noticing that this approach gives the same values of the crit- 
ical temperature for the two- and three-dimensional Ising ferromagnets as the lowest 
linear approximation within the modified effective-field theory [13] in the so called P 
scheme. Consequently, the present EFT yields, for example, a nonzero critical con- 
centration for quenched diluted systems, the lack of order in the one- dimensional 
Ising ferromagnet, and the occurrence of order in the two-dimensional case with 
the critical temperature improved over the usual mean-field theory [13] (for another 
advanced effective- field methods see also, e.g., Ref. [15]). For the present frustrated 
pure system, in contrast to the mean-field theory, our simple EFT approach [TO] 
correctly reproduced no long-range order behavior down to zero temperature in zero 
field and a fairly accurate phase diagram in finite fieldfl 



[Fig. 1 about here.] 

In the present paper, we generalize the study done in Ref. [TU] by consider- 
ing the concentration p sublattice-dependent and hence characterized by a three- 
component vector p = (pa,Pb,Pc)- Following the procedure described in Ref. [10J, 
the A-sublattice magnetization per site in zero field can be calculated by taking a 



Comparison can be made with some other more accurate techniques, such as Monte 
Carlo and renormalization-group studies [3H5] 
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configurational average of the expression 



JAB 



^(Sh =C A ( II cosh (/3 JD) + $Sf sinh(/3 JD) + 1 - (2) 
x n [4° cosh(/3 JD) + 4 c 5 fc c sinh(/3 J/J) + 1 - ) tanh(x) \ x=0 , 

k=l ' 

where z AB ,z A c are the numbers of NNs of the spin from the sublattice A that 
belong to the sublattices B and C, respectively, /3 = 1/ksT and D = d/dx is the 
differential operator. Similar expressions for the sublattices B and C can be obtained 
from Eq. (j2J) by cyclic permutation of the indices A, B and C and for each sublat- 
tice considering an appropriate number of NNs on the remaining sublattices (in our 
case zxy = 3, for all pairs X, Y = A, B or C). Then, by the use of the decoupling 
approximation for the thermal averaging [TU] and performing the configurational av- 
eraging of these equations for the respective sublattices leads to a system of coupled 
equations for the averaged sublattice magnetizations per site in the form 

m A = Pa (ob + bm B f (a c + bm c f tanh (x) \ x=0 , 

m B = p B (ac + bmcf (a A + bm A ) 3 tanh (x) \ x=0 , ( 3 ) 

m c = Pc (a>A + bm A f (a B + bm B f tanh (x) \ x=0 , 

where ax = 1 — Px + Px cosh (f3JD), with X = A, B or C, and b = smh(flJD). 
Using the differential operator relation exp(aD)f(x) = f(x + a), Eqs. (j3J) can be 
expressed explicitly 

In zero field, the critical temperature can be determined as follows. Using the fact 
that in the vicinity of the second-order phase transition the sublattice magnetizations 
mx (X = A, B or C) are very small, Eqs. ([3]) can be expanded and linearized, which 
leads to a system of homogeneous linear equations 

Um = 0, (4) 

where 



— 1 3p A a B a c b tanh (x) \ x=0 3p A a B a c b tanh (x^ 1 ^ 



U 



\x=0 
\x=0 



3pBa A a c b tanh. (x) \ x= o —1 3pB«i a c^tanh (x 

3pca\a B btanh(x) \ x= o 3pca A a B btanh(x) \ x= o —1 



(5) 



and m = (mA,mB,m c ) T . The critical (or Neel) temperature T N can be established 
as a function of the sublattice concentration vector p = (p A ,PB,Pc) by finding a 
nontrivial solution of Eqs. (jl]), i.e., by solving the secular equation 

detU = 0. (6) 
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3 Results and discussion 

In the following, we will consider a selectively diluted case such that two (let us say 
A and B) sublattices are diluted while the remaining sublattice (C) is occupied with 
purely magnetic ions, i.e., the case of p = (j>a,Pb,1), where < Pa,Pb < 1- The 
special case of p = (pa, 1, 1), when only one sublattice is diluted, was investigated by 
hard-spin mean-field theory (HSMFT) [11] and Monte Carlo (MC) simulations [12J. 
Compared to the regular mean-field theory, HSMFT conserves frustration since, 
instead of a uniform magnetization, it considers local magnetizations at each of N 
sites which feel the antialigning field due to the full (i.e., hard) spin of each of its 
neighbors. A further approximation (FA) of HSMFT still incorporates frustration, 
however, the approach is simplified by imposing sublatticewise uniformity. Thus, FA 
HSMFT leads to a set of only three instead of N coupled equations. These HSMFT 
and MC simulation works have shown that below some threshold concentration 
p% the nondiluted sublattices develop a long-range order (LRO) with nonvanishing 
magnetizations m-Q = —uiq that do not saturate as the temperature approaches 
zero. On the other hand, the magnetization of the diluted sublattice tua remains 
zero at any concentration. 

[Fig. 2 about here.] 
[Fig. 3 about here.] 

A similar picture is also observed in our EFT approach. The resulting phase diagram 
is plotted in Fig. [21 in which the area below the curve represents the LRO phase. This 
diagram can be compared with similar plots obtained from HSMFT and FA HSMFT, 
shown in Fig. 3 of Ref. [UJ. In both EFT and HSMFT approaches the transition 
temperature decreases with increasing concentration, albeit, the EFT results seem 
to be quantitatively closer to those obtained from FA HSMFT rather than HSMFT. 
Namely, the EFT critical value of p c A = 0.8753, above which no LRO exists in either 
of the sublattices, matches well with the one obtained from FA HSMFT (p A = 0.875) 
but it is lower than those obtained from HSMFT (p c A = 0.958) and MC simulations 
(p c A ~ 0.95) [12], respectively. In the limit of pa = the EFT critical tempera- 
ture of the honeycomb lattice model of feeT^/l J\ = 2.1037 is accurately recovered 
(see, e.g., Ref. JTH] where also comparison with some other methods is provided). In 
Fig. [3] we plot a temperature variation of the nondiluted sublattice magnetizations 
ins = —mc for different values of p\. By comparing it with Fig. 1(a) of Ref. [TTj . 
one can see a qualitatively similar behavior: the sublattice magnetizations do not 
saturate at zero-temperature, except for the full dilution (honeycomb lattice) limit, 
and not only the transition temperatures (as already seen in Fig. [2] shown above) 
but also the values of the sublattice magnetizations at zero temperature decrease 
with increasing concentration pa- 

Now let us see how this picture will change if the dilution is extended to two 
sublattices, i.e., if p = (pa,Pb, 1), where < Pa,Pb < 1- The corresponding phase 
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diagram in the (pa, Pb) plane is presented in Fig. HI We note that with respect to the 
sublattice magnetizations, the critical surface shown in Fig. H] represents the tem- 
peratures above which all the sublattice magnetizations vanish. However, in case 
of p = (pa, 1, 1) presented above, the system can display a partially ordered phase 
with one (diluted) sublattice magnetization equal to zero even below the critical 
surface. The above mentioned special case of p = (pa, 1,1), shown in Fig. El is 
represented by the curve of pb — 1 and varying pa- When we also start diluting 
the sublattice B we observe overall decrease of the transition temperature and until 
Pb = Pa also decrease of the value of p c A from 0.8753 at p B = 1 down to 0.7640 
at pb = Pa- Moreover, we have found that for p B 6 (0.7640, 0.8753) the system 
shows the reentrant behavior. Namely, at a fixed value of pb by varying the value of 
Pa the system passes from the ordered phase to the disordered one and then again 
to the ordered one. This behavior can be more clearly seen in Fig. [5] in which the 
ground-state phase diagram is depicted. The latter can be obtain from Eq. ([6]) in 
zero-temperature limit. For p B G (0.4280, 0.7640), LRO persists for any value of pa, 
however, at pb = Pa the transition temperature curve changes from decreasing to 
increasing. For pb £ (0.2976, 0.4280) the system again displays the reentrant behav- 
ior with LRO disappearing within a certain range of the concentration pa- Upon 
further decrease of pb there is no LRO at lower values of pa but it appears at some 
threshold value p A , which decreases down to p c A = 0.2976 at pb = 0. The limiting 
case of pb = (or equivalently pa = 0) corresponds to the honeycomb lattice. There- 
fore, if we consider the case of pa = 0, the dependence of the critical temperature 
on pb represents the phase diagram of the selectively diluted Ising antiferromagnet 
on the honeycomb lattice and the value p B = 0.2976 is the percolation threshold, 
i.e., the critical concentration below which LRO cannot exist due to the absence 
of a spanning cluster of magnetic ions (Fig. [6]). We note that the phase diagram is 
symmetric along the diagonal and therefore the above arguments also apply when 
Pa and pb are interchanged. 

At this point, it is interesting to examine the frustration effect on the perco- 
lation behavior by comparing it with its nonfrustrated counterpart, i.e., with the 
ferromagnetic case. The selective site percolation problem in an Ising ferromagnet 
on triangular and square lattices with several kinds of sublattice structures has been 
extensively studied by Idogaki and Uryu [T714T9] , employing a renormalization group 
technique. However, it is difficult to assess the frustration effect by a direct compari- 
son of the percolation thresholds from our study with those obtained in Refs. [TTlfrQ] 
due to generally different sublattice structures considered as well as different tech- 
niques used. Therefore, we made additional calculations for the current selectively 
diluted model within EFT but considering the exchange interaction J > 0, i.e., for 
the triangular Ising ferromagnet, and the results are presented by the dashed curve 
in Fig. [5J Naturally, for the ferromagnetic case there is no disordered region at high 
concentrations and the low concentration disordered region is much smaller than 
for the frustrated antiferromagnetic case. As can be expected, the biggest difference 
between the two cases is along the diagonal (j>a = Pb), where the frustration in the 
former case is the highest, and vanishes in the limits of pa = and pb = 0, where 
there is no frustration. 
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[Fig. 4 about here.] 
[Fig. 5 about here.] 
[Fig. 6 about here.] 

For the case of p = (p\, 1, 1), investigated in Refs. [TT|[T2] as well as in the current 
study, it has been shown that below the threshold value of p\ the nondiluted sublat- 
tice magnetizations become finite and the diluted sublattice magnetization always 
remains zero. Moreover, at zero temperature the nondiluted sublattice magnetiza- 
tions do not saturate for any nonzero concentration. In order to shed more light 
on these phenomena, let us study the ground-state behavior of the current model 
with two sublattices diluted, i.e., the case of p = (pa,Pb, 1)- It is well known that in 
the pure model the elementary triangular plaquette in the ground state is six-fold 
degenerate. The selective site dilution lifts the degeneracy and the system can show 
long-range order. The possible zero-temperature sublattice arrangements and their 
respective energies are presented in Table HJ It is easy to verify that if both sublat- 
tices are diluted, i.e., Pa,Pb < 1, then the lowest energy is e 3 = PaPb — Pa — Pb, 
and the corresponding configurations are those with spins on the diluted sublat- 
tices A and B oriented in one direction and spins on the nondiluted sublattice C 
in the opposite direction. On the other hand, if only one sublattice is diluted, e.g., 
Pa < 1 and pb = 1, then the lowest energy is e± = = — 1, and the corresponding 
configurations are those with spins on the nondiluted sublattices B and C oriented 
antiparallel but spins on the diluted sublattice A can point in either direction with 
equal probability. This four-fold degeneracy explains why the diluted sublattice mag- 
netization remains zero. 



[Table 1 about here.] 

The unsaturated magnetizations of the nondiluted sublattices for the p = (pa, 1,1) 
case, shown in Fig. [3j apparently result from the presence of the frustration. As the 
dilution (frustration) increases (decreases), the magnetizations gradually reach the 
saturation value of 1. The initial increase is fast and for pa = 0.5 about 96% of 
the saturation value is reached. Nevertheless, with decreasing pa the curve flattens 
out and the full saturation is only achieved at pa = 0. This behavior is illustrated 
in Fig. in which we plot the sublattice magnetizations at zero temperature. The 
results are in good agreement with those presented in Fig. 2(a) of Ref. [UJ. Further, 
due to considerable variation of the frustration degree in the (pa,Pb) parameter 
space, it will be interesting to see similar plots also for the p = (pa,Pb, 1) case at 
different values of pa and p-Q- The reached saturation level is better seen if instead 
of the sublattice magnetizations per site nix we plot the sublattice magnetizations 
per spin fix = nix/px, where X = A,B or C. Note that in Fig. [7J fix = m>x- 111 
Fig. [8] we present the variations of the zero-temperature sublattice magnetizations 
per spin in different regions of the (pa,Pb) plane. For p B = 0.8, as expected from 
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Figs. H] and [5], we observe the reentrant behavior featuring a certain region of the 



values of pa in which /i x = 0, for X = A,B or C (Fig. 8(a) ). Outside of this highly 
frustrated region, in accordance with the ground-state analysis presented in Table [TJ 
the diluted sublattices A and B have nonzero magnetizations with the signs opposite 
to that of the nondiluted sublattice C, and in the limit of pa = 1 the magnetization 
of the remaining diluted sublattice B reduces to zero. Interestingly, in the other 
limit of pa — >■ 0, although the A sublattice magnetization per site rriA vanishes, the 
value of ha remains finite and drops discontinuously to zero only at pa = 0. Also, 
the magnetization of the nondiluted sublattice C falls short of saturation. A closer 



inspection of Fig. 8(a) reveals that the sublattice magnetization /ic does not reach 
its full saturation value of -1 even at pa = 0, i.e., for the selectively diluted hon- 
eycomb lattice system. Thus, the lack of saturation in this case cannot be ascribed 
to the frustration but rather to the presence of small amount of isolated C-ions or 
isolated clusters involving both C- and B-ions surrounded by vacancies which do 
not contribute to the magnetizations of the C- or both C- and B-sublattices, respec- 
tively. The above described phenomena can be observed even more pronouncedly 



for pb = 0.4 (Fig. 8(b)). The more abundant presence of the "loose" clusters in- 



volving also B-ions with decreasing pb is reflected in the decreasing value of /iB at 



Pa = 0. Here we note that in both Fig. 8(a) and (b) the source of the unsaturation 



in the limit of pa = is different from that observed in the limit of pa = 1- In the 
former case the system is nonfrustrated (selectively diluted honeycomb lattice) and 
the lack of saturation can be ascribed solely to the presence of the "loose" isolated 
C-ions or clusters of C- and B-ions, which do not align with the spanning clusters. 
On the other hand, in the latter case the system is frustrated (selectively diluted 
triangular lattice) and the unsaturation appears in the nondiluted sublattices. Since 
the random dilution of the A-sublattice relieves frustration locally, the alignment of 
local magnetizations in the nonfrustrated sublattices proceeds in a gradual manner. 
Therefore, one can expect the presence of some amount of antialigned local magneti- 
zations that cause the lack of saturation until the A-sublattice is completely removed 
by dilution. For the intermediate concentrations < pa < 1 the unsaturation effect 
is amplified by the presence of both sources. Namely, by increasing (decreasing) pa 
from zero (one) the increasing frustration (dilution) magnifies the degree of unsat- 
uration until the long-range order is destroyed in all the sublattices. In Fig. 8(c) we 
show the case of pb = 0, i.e., the selectively diluted honeycomb lattice case. In this 
case the unsaturated magnetizations are due to isolated C-ions and clusters of C- 
and A-ions. If the sublattices A and B are diluted with equal probabilities, then the 
frustration is the least relieved and LRO with considerably suppressed values of the 
zero-temperature sublattice magnetizations appears only in a relatively small range 



of 0.4280 < pa = Pb < 0.7640 (Fig. 8(d)). On the other hand, the frustration is the 



most avoided if the dilution is performed such a way that pa + Pb = 1 (Fig. 8(e)) 
In this case LRO survives at any dilution. 

[Fig. 7 about here.] 

[Fig. 8 about here.] 
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4 Summary and conclusions 

We studied the critical and ground-state behavior of a selectively site-diluted Ising 
antiferromagnet on triangular and honeycomb lattices within the framework of an 
effective-field theory with correlations. The selective dilution was carried out by ran- 
dom removal of magnetic ions from two of the three sublattices with different prob- 
abilities. In a special case when only one sublattice is diluted, our results showed a 
fairly good agreement with the previous studies by hard-spin mean-field theory [11] 
and Monte Carlo simulations [12J. The extension to the system with two sublat- 
tices jointly diluted with different probabilities is an entirely new contribution. We 
established a three-dimensional phase diagram in the space of the sublattice concen- 
trations of the triangular lattice. We identified one smaller disordered region at high 
sublattice concentrations in which a long-range order cannot occur due to high frus- 
tration. At low sublattice concentrations there is another disordered region in which 
a long-range order is absent due to the absence of a spanning cluster of magnetic ions. 
The two disordered regions are stretched along the diagonal owing to the increased 
frustration when the sublattice concentrations are about the same. These relatively 
small regions are surrounded by one larger region in which the sublattice concen- 
trations are large enough to form a spanning cluster and the frustration is relieved 
enough to enable long-range ordering. We also studied the ground-state ordering and 
the effect of unsaturated sublattice magnetizations. The degree of unsaturation was 
found to be proportional to both the degree of frustration and dilution. However, 
due to compensating effects of the two factors (the increasing dilution decreases frus- 
tration), the highest unsaturation was observed at the intermediate dilution range, 
in the vicinity of the phase boundaries. The selectively diluted Ising antiferromagnet 
on the honeycomb lattice was obtained as a special case when one sublattice of the 
triangular lattice is completely diluted. We established the phase diagram including 
the percolation threshold below which no long-range order can exist. 

In the present study we demonstrated how the selective dilution of one or two 
sublattices in the triangular lattice antiferromagnet can relieve frustration and cause 
long-range ordering, featuring some peculiar phenomena, such as reentrance and un- 
saturated sublattice magnetizations. An introduction of an external magnetic field 
is another way of perturbation of the highly degenerated triangular lattice antiferro- 
magnetic system and, therefore, we believe that it would be interesting to perform a 
similar study of the combined effects of the selective dilution and the external field. 
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Sublattices: A-O B - □ C-A 

Fig. 1. Triangular lattice partition into three sublattices A,B and C. The shaded sym- 
bols signify the selective dilution case when two sublattices A and B are diluted, i.e., 
P = (PAiPB, !)• 
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k B T/|J| 

Fig. 3. Temperature variation of the nondiluted sublattice magnetizations m-Q = -mc, for 
different values of pa, for the case p = (p\, 1,1). 
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Fig. 4. Transition temperature as a function of the concentrations px and pb, for 
P = (PAiPB, 1)- The bold lines outline the zero-temperature limit. 
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Fig. 5. Ground-state phase diagram in p\ — pb concentration plane, for p = (pa,Pb, !)• 
LRO and D denote the long-range ordered and disordered phases, respectively. The dashed 
curve corresponds to the nonfrustrated ferromagnetic case. 
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Fig. 6. Transition temperature as a function of the concentration p-Q, for p = (0,pb,1)- 
This dependence represents the phase diagram of the selectively diluted antiferromagnet 
on a honeycomb lattice. 
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Fig. 7. Zero-temperature sublattice magnetizations for p = (p\, 1, 1). 
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Fig. 8. Zero-temperature sublattice magnetizations per spin for (a) pb = 0.8, (b) pb = 0.4, 
(c) pb = 0, (d) pa = PB, and (e) p A = 1 - p B . 
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Table 1 

Zero-temperature degeneracy with possible sublattice arrangements (tua, mB, m c) an d 
corresponding energies per site (H)/\J\N for p = (pa,Pb, !)• 



n 


12 3 


(mA,m Bl mc) 


(+,-,+);(-,+,-) (-,+,+);(+,-,-) (-,-,+);(+,+,-) 


(H) 

&n ~ \J\N 


-PAPB +PA-PB -PAPB -PA+PB PAPB ~ PA ~ PB 
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